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ABSTRACT 

We model the multi-wavelength emission in the southern hotspot of the radio quasar 4C74.26. 
The synchrotron radio emission is resolved near the shock with the MERLIN radio-interferometer, 
and the rapid decay of this emission behind the shock is interpreted as the decay of the amplified 
downstream magnetic field as expected for small scale turbulence. Electrons are accelerated to 
only 0.3 TeV, consistent with a diffusion coefficient many orders of magnitude greater than in 
the Bohm regime. If the same diffusion coefficient applies to the protons, their maximum energy 
is only ~ 100 TeV. 

Subject headings: galaxies: active — galaxies: jets — quasars: individual(4C74.26) — acceleration of 
particles — radiation mechanisms: non-thermal — shock waves 


1. Introduction 


Diffusive shock acceleration (DSA) is an estab¬ 
lished mechanism to convert bulk kinetic energy 
into a non-thermal distribution of relativistic par¬ 
ticles with a maximum energy much larger than 
the average energy of particles in the plasma. This 
theory explains well the spectrum of Galactic cos¬ 
mic rays (CR) with energies up to PeV, ac¬ 
celerated in supernova remnant shocks (see 


Bell 


2014 for a review). The most energetic CR, i.e. 


particles with energies up to 50 EeV, are acceler¬ 
ated outside the Galaxy but the origin of these 
particles is still unknown. Relativistic shocks 
in extragalactic sources like Gamma Ray Bursts 
and Active Galactic Nuclei ha ve been pro posed 
as candidates (e.g . iGallant fe Achterberd Il999l 


Murase et al.ll2012f ). However^theoretical models 


( Pelletier^taL 20091 Lemoine fe Pelletier 2010t 
Sironi et al.1 2013l ~ Reville fe Bell 2014 ) show mag¬ 


netic field amplification at ultra-relativistic shocks 
on scales much smaller than the Larmor radius 
r g of particles being accelerated which precludes 
CR acceleration to EeV energies unless other pro¬ 
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cesses can be found to amplify the magnetic field 
on larger scales. 

Hotspots are usually detected at the jet termi¬ 
nation re gion in type II Fanaroff- Riley (FR) radio¬ 
galaxies ( Fanaroff fc Rilev 1974 ). The location of 
the hotspot is coincident with the downstream re¬ 
gion of the jet termination shock, where particles 
accelerated by the shock emit synchrotron radi¬ 
ation. Therefore, hotspots are suitable places to 
study DSA in high velocity shocks. 

We model the emission from radio to X-rays in 
the southern hotspot of the FR II source 4C 74.26 
using data provided in Erlund et al. ( 201Clh . We 
determine that the compact radio emission traces 
out the location of the shock where the magnetic 
field B is amplified by plasma instabilities up to 
~100 /iG, and it damps rapidly downstream of the 
shock. The turnover in the synchrotron spectrum 
between infrared (IR) and optical wavelengths im¬ 
plies that the CR electron scattering length is 
much longer than the Larmor radius and consis¬ 
tent with the amplified magnetic field being struc¬ 
tured on very small scales and comparable with 
the ion skin-depth c/w p ;. Cosmic ray acceleration 
is consequently very slow and electrons are accel¬ 
erated to only E Stmax ~ 0.3 TeV. A similarly low 
acceleration rate for ions would limit their energy 
to ~ 100 TeV and preclude proton acceleration to 
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EeV energies. 


2. The giant FR II galaxy 4C74.26 


The FR II galaxy 4C74.26 is located at red- 
shift z = 0.104 (^0.5 Gpc from EarthJ^. Two 
X-ray sources sep arated by ~10" we re detected 
with Chandra (by lErlund et al. 201 0 : ) at the ter¬ 
mination region of the southern jet, as shown 
in Figure [Q (upper). In the present work we 
study the southern X-ray source, with a luminosity 
L x ~ 10 41 erg s _1 at 2 keV and called “the south¬ 
ern arc”. The shape of this emission is arc-like 
with a characteristic size l x ~ 10", and encloses a 
compact radio source. 


2.1. The southern arc 


Compact radio emission from the southern 
arc was detected with the MERLIN high reso¬ 
lution interferometer (z/ r = 1.66 GHz) with a flux 
/ r ~0.04 Jy, and luminosity L r ~ 1.9xl0 40 erg s _1 
per unit logarithmic bandwidth 5v = v. This 
emission is located in a region of width Z r < 1" on 
the plane of the sky. In addition, faint and diffuse 
radiation was detected at IR (i/; r = 1.36xl0 14 Hz) 
and optical (i/ opt = 6.3 xlO 14 Hz) bands, with 
fluxes ~8.4xl0 -6 and 2.82xl0~ 7 Jy, respectively, 
and located in a region of width > l r . However, 
there is a linear structure (in both bands) that 
traces the brightest edge of the MERLIN radio 
emission, and seems to be cupped within it. 


Two factors indicate that the southern arc of 
X-ray emission is not synchrotron. First, l x > l r 
is inconsistent with the X-ray emitting electrons 
being more energetic and therefore cooling more 
rapidly as they advect away from the shock. Sec¬ 
ond, the steep spectrum between IR and optical 
(see Fig. 13 in lErlund et, al.1 l2010h indicates the 
maximum energy of (synchrotron) emitting elec¬ 
trons. We note that similar characteristics are o b- 
served in other sources (e.g. lOrienti et al. 2012 1. 

lErlund et al. ( 2010h suggested that the multi¬ 
wavelength emission from the southern arc is pro¬ 
duced by non-thermal electrons, emitting syn¬ 
chrotron radiation from radio to optical, and 
up-scattering the cosmic microwave background 


1 Throughout this paper we use the cosmology Ho = 
71 km s —1 Mpc~', fl 0 = 1 and Ao = 0.73. One arcsecond 
represents!.887 kpc on the plane of the sky at z = 0.104. 




Fig. 1.— Upper: the hotspot(s) at the southern 
jet of the FR II gala xy 4C 74.26 (a dapte d and ro¬ 
tated by ~60° from lErlund et al.l 12010 1. White 
and yellow contours are X-rays and radio data, 
respectively. Red and green correspond to IR 
and optical, respectively. Lower: sketch of our 
model of the southern arc (not to scale). The 
synchrotron radio-to-optical radiation is located 
within the compact MERLIN emitter, whereas IC 
X-ray emission is produced in a more extended 
region. 
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(CMB) photons (with energy ~7xl0 -4 eV and 
energy density U C mb = 6 xl 0 -13 erg cm -3 ) to the 
X-ray domain. In this scenario, the radio-to-IR 
spectral index is a = 0.75 (the synchrotron flux 
density at frequency v is /„ oc v~ a \ which cor¬ 
responds to p = 2.5 in the non-thermal electron 
energy distribution when the emission is produced 
in the same volume. 


3. The hotspot as a magnetic field damp¬ 
ing region 

In this work we con sid er t h e sa me emission 
mechanisms as in lErlund et al. ( 2010h . but allow 
the synchrotron and Inverse Compton (IC) emis¬ 
sion to be produced in regions with different spa¬ 
tial extents. In particular; electrons accelerated at 
the shock emit synchrotron radiation from radio 
to optical in a compact region behind the shock, 
whereas the IC X-ray emission is located in an 
extended region. In Fig. |T] (lower) we sketch our 
model. 

The synchrotron (s) and IC cooling length 
of electrons with Lorentz factor 7 is 4,ic(b) = 
t s ,ic( bVsh/c where t s> i c (b) is the cooling timescale. 
The shock velocity is approximately the same as 
the jet velocity which we take characteristically 
to be u s h = 10 10 cm s -1 (~ c/3 and Lorentz fac¬ 
tor r sh ~ 1.06) in line with observations o f similar 
objects (see lSteenbrugge fe Blundell 2008 . and ref¬ 


erences therein). We use r = 7 as the shock com¬ 
pression ratio for a non-relativistic shock whose 
downstream thermal pressure is dominated by rel¬ 
ativistic electrons, although r ~ 4 may still ap¬ 
ply if non-relativistic ions domin ate t h e pr essure 
downstream of the shock ( Kirk et al. 200 0). Our 
conclusions are not sensitive to the exact value of 


3.1. Inverse Compton X-ray emission 

The IC X-ray emission is produced by electrons 
withyx ~ 10 3 and ?i c (7x) ~ 10 4 (v s h/10 10 cms -1 ) arc- 
sec, which is much larger than Z x . The synchrotron 
cooling length Z S ( 7 X ) is also greater than Z x , unless 
the magnetic field in the X-ray emitting region is 
~ 360 /iG. However, such a large magnetic field 
would produce synchrotron radio radio emission 
much brighter than L r (see next section). Further¬ 
more, we show below that the amplified magnetic 
field, of the order of ~100 /iG, is confined to a 


small volume close to the shock. Therefore, adia¬ 
batic expansion is probably the dominant cooling 
mechanism as the particles flow out of the hotspot. 


Unless F s h > 10, X-ray emitting electrons are 
non-thermal and fo llow a p ower law ener gy distri¬ 
bution oc 7 -p (e.g. Giannios fe Soitkovskv 2009h . 
Assuming that the X-ray emitting volume is 14 ~ 
300 arcsec 3 , the energy density of these non- 
thermal electrons is ^10 _9 (7 m i n /50) _o - 5 erg cm -3 
where the power law terminates at a minimum 
Lorentz factor 7 m i n . The magnetic field with the 
same energy density is ~ 100(7 m i n /50) -0 - 25 /tG. 
These results correspond to the case where p = 2.5 
(see however Sect. 13.2221) . In the following sections 
we take 100 pG as a fiducial magnetic field since 
it represents equipartition between magnetic and 
relativistic electron energ y densities and is typic al 
of other hotspots (e.g. Godfrey fc Shabala 2013 b 


3.2. Synchrotron emission 

Considering that 7 ( 1 /) ~ 4.5 x 10 
is the Lorentz factor of electrons emitting syn¬ 
chrotron radiation at v in a magnetic field B , l B 
can be written as 

l Y) _ 12 ( B A" 1 ' 5 ( v sh A 

[//] VGHz/ \100/tG J \10 lo cms -1 / 

(1) 


3.2.1. Radio 

MERLIN emitting electrons have 7 r = 7 ( 27 ) ~ 

27 x (H /100 /iG) -0 - 5 . If both radio and X-ray emis¬ 
sion are produced by non-thermal electrons that 
follow the same power-law energy distribution, 

LJL r ~ (7x/7r) 3 “ P (^cmb/17mag)K/f 4 r and B ~ 

100 /iG corresponds to 14 /V T ~ 5 x 10 3 , where 
U x nag = B 2 /8n and V T is the volume of the 
synchrotron emittei0. Such a large ratio be¬ 
tween emitting volumes is not implausible pro¬ 
vided the magnetic field is inhomogeneous in the 
shock downstream region and the synchrotron 
emitter consists of features smaller than the MER¬ 
LIN point spread function (FWHM 0.15") as seen 
in parts of the MERLIN data. 

The synchrotron cooling length of MERLIN 
emitting electrons is l a (v r ) ~ 9.3"(.B/100 /iG) -3 / 2 (i; s h/10 10 cms -1 ), 
and a very large magnetic field of ~ 2.4 (u s h/10 10 cm s -1 ) 2 / 3 mG 


2 Note that if V x = V r , an unrealistically small magnetic field 
of 0.3 nG would be needed to explain the observed fluxes. 
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would be required to match Z a (i/ r ) = 0.1" ~ l T . 
This result suggests that the downstream extent 
of the compact emission detected at v r is not the 
result of fast synchrotron cooling, as we can con¬ 
firm when we take into account the IR emission. 

3.2.2. Infrared 

The synchrotron cooling length of IR emitting 
electrons is Z s (^ir) ~ 0.03"(R/100/xG) _3 / 2 (u s h/10 10 
indicating that these particles radiate most of their 
energy within l T . (This angular distance is not re¬ 
solved by the IR observations.) This is consistent 
with a radio-to-IR electron energy spectral index 
of p ~ 2.5 with the cooling break in the spectrum 
occuring close to IR wavelength^. Note that if the 
emitting volume were determined by synchrotron 
cooling, l s (u) oc v~ 0 5 giving p = 2a = 1.5 since a 
is measured to be 0.75. This very hard spectrum 
is unlikely since it diverges toward high energy 
and would be remarkable in hotspots, supporting 
the conclusion that the downstream radio extent 
l T must be determined by factors other than syn¬ 
chrotron cooling. As we discuss in Sect. |4] this 
may be t he resu lt of the dam ping of the magnetic 
field (see lSchure et al.|[2012i for a review). 

3.2.3. Optical 

Optical emission produced by synchrotron ra¬ 
diation of electrons with 7(r , op t) ~ b(r'ir) is al¬ 
most co-spatial with the IR emission, and this 
explains the linear structure cupped within 
The synchrotron turnover v c between v i r and 
i/ op t indicates that the maximum energy of non- 
thermal electrons is E e ,ma.x ~ 'y{i / c )m e c 2 ~ 
0.3(i/ c /^ r ) o - 5 (B/100AtG)-’°- 5 TeV. 


was provided by iBell ( 20041) showing that non¬ 
resonant hybrid instabilities are capable of en¬ 
hancing the magnetic field by orders of magnitude. 
Magnetic field amplification is also responsible 
for B ~ 100 pG in the southern arc in 4C74.26 
since it is much larger than the expected value 
in the jet u pstream o f the termination shock (e.g. 
Hardcastle fc Krausell2014l ). 


cms 


-i -j Bohrn diffusion (electron mean free path A ~ 
9 - g ) in a ~100 pG magnetic field would be ex¬ 
pected to accelerate electrons with synchrotron X- 
ray emittin g energies as seen in supernova rem¬ 
nants ( Stage et all 2006). However, v c ~ r'ir.opt 
determined by a competition between shock ac¬ 
celeration and synchrotron cooling indicates that 
acceleration is slow and therefore that the elec¬ 
tron diffusion coefficient D is much larger than 
the Bohm value Dsohm: 


D 


10 b 


V s h 


O’G 


( 2 ) 


A>Bohm V10 10 cms 

independent of B (see e.g. Casse et ~ad]l2002h . Such 
a large diffusion coefficient in an amplified mag¬ 
netic field is expected if it is structured on a scale 
s much smaller than the Larmor radius of the elec¬ 
trons being accelerated. Small angle scattering by 
magnetic field randomly orientated in cells of size 
s produces D 


— ~ 10 7 ( — 

cm 


(r g /s)£>Bohm and then 

(: 


g 

1.5 


B 


100 /zG 


- 1.5 


^sh 


. 10 10 cms -1 

( 3 ) 

In comparison the ion skin-depth is c/w p ; ~ 2.3 x 
10 9 (n/10 -4 cm -3 )” 0 - 5 cm, where n is the particle 
density downstream of the shock (assumed to be 
7 times the jet density), and 


4. Magnetic field amplification 


The amplification of the magnetic field at 
strong sh ocks in supernova remn ants was demon¬ 


strate d bv IVink fc Laming ( 2003 ) and Berezhko et al 


( 20031 ). deriving the magnetic field from l s down¬ 
stream of the shock. A theoretical explanation 


3 Note that the relationship p = (r + 2)/(r — 1) b r eaks down 
for m ildly relativistic shocks iKirk et all 2000; Bell ct al. 
20li|) or when non -linear feedback is important (e.g. 
Amato fc Blasill2005l ). 

4 The faint diffuse IR and optical emission may be the result 
of CMB photons up-scattered by electrons with 7 ~ 50. 


s 

O.O! 

y. 5 

' «sh \ 2 

c/CJpi 

) 1 

k 10 1U cms -1 ) 


( B ) 

-1.5 

( n l 0 - 5 


V 100 mG J 


V 10 -4 cm -3 J 


Given the uncertainties in the parameter val¬ 
ues, the approximate nature of the theoreti¬ 
cal models, and the wide range of the spatial 
scales (s, r g , Z r ), it is not significant or surpris¬ 
ing that our estimate of s/(c/w p i) differs from 
unity by a factor of ~ 0.01. The order of magni¬ 
tude similarity of s and c/oj p , supports the con¬ 
tention that shock-generated small-scale turbu¬ 
lence scatters non-thermal electrons during diffu¬ 
sive shock acceleration. This is consistent with 
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Sironi &_Spitkovsky (2011) who discuss the var¬ 
ious processes related to the Weibel instability 
that excite turbulence on the characteristic scale 
of c/ujpi. Simulations show that magnetic field 
generated by the Weibel instability decays down¬ 
stream of the shock because of it s rela tively small 
scale l ength (Sironi fc Spitkovskvll201ll : I Bret et aT 


2013; Sironi et ah 2013 ). This would account for 
the cut-off of synchrotron emission in 4C74.26 far 
short of the synchrotron cooling distance of radio- 
emitting electrons (l T <C l s (y r )). These electrons 
continue up-scattering CMB photons, thus pro¬ 
ducing IC X-ray emission downstream of the shock 
after the MERLIN radio emission has ceased. 


4.1. Limit on ions maximum energy 

Since the response of highly relativistic ions is 
similar to that of electrons with the same energy in 
a tangled amplified magnetic field, we can expect 
protons to have a similar ratio of -D/Dsohm- Pro¬ 
tons can be accelerated to higher energies than 
electrons because their radiative losses are min¬ 
imal, but the maximum energy to which they 
are accelerated is limited because their acceler¬ 
ation time is increased by t he rati o -D/Dsohm- 
The Hillas parameter v B h B R ( Hillas 1984h . where 
R ~ 2" (~ 2.8 kpc) is the characteristic length 
of the source, would suggest proton accelera¬ 
tion to ~100 EeV in the termination shock of 
4C74.26, but the maximum energy is reduced to 
only ~100 TeV if D ~ 10 6 Dsohm since the Hillas 
parameter assumes D ~ HBohm and is otherwise 
reduced by the factor (H/HBohm) -1 - Another 
perspective on the same effect is that the mean 
free path for scattering by small-scale turbulence 
A ~ r 2 g /s is larger than the size of the system if 
s ~ c/wpi and r g is the Larmor radius of an EeV 
proton. This result suggests that the mildly rel¬ 
ativistic termination shock in 4C74.26 is a poor 
accelerator of UHECR. 


be determined by synchrotron cooling through this 
wavelength range (Sect. [3T2~2ll . 3) The turnover of 
the synchrotron spectrum at IR/optical frequen¬ 
cies requires A>r g for any reasonable shock ve¬ 
locity (Sect. HI). These three features fit well in 
a scenario in which the MERLIN radio emission 
traces out the region where the magnetic field is 
amplified by plasma instabilities with small length 
scale (e.g. Weibel). 

The magnetic field in equipartition with non- 
thermal electrons in the MERLIN emission re¬ 
gion is ~100 uG and simila r to t he val ues ob¬ 
tain ed by other authors (e.g. Godfrey & Shabala 
20131 ). An unrealistically large magnetic field 


2.4 (u s h/10 10 cms -1 ) 2 / 3 mG would be needed 
to explain the compact radio emission in terms of 
synchrotron cooling. If B ~ 100 /iG in the syn¬ 
chrotron emission region, the maximum energy of 
non-thermal electrons is ~0.3 TeV. 


If ions are accelerated as well, protons with en¬ 
ergy ~0.3 TeV diffuse also with mean free path 
A > r g . If A is similarly larger than the Larmor 
radius at higher proton energies, then the maxi¬ 
mum proton energy at the termination shock of 
4C74.26 is only 100 TeV instead of the 100 EeV 
indicated by the Hillas parameter. This may have 
important implications for the understanding of 
the origins of UHECR. 
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